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Abstract. We present the results of 20 ksec Chandra observations for each of 5 radio galaxies in the redshift range 
2.0 < z < 2.6. The goals were to (i) study the nature of their non-thermal X-ray emission, (ii) investigate the 
presence and amount of hot gas, and (iii) look for active galactic nuclei (AGN) overdensities in helds around high 
redshift radio galaxies. For 4 of the 5 targets we detect unresolved X-ray components coincident with the radio 
nuclei. From spectral analysis of one of the cores and comparison to the empirical radio to X-ray luminosity ratio 
(Lr/Lx) correlation for AGN, we find that the cores are underluminous in the X-rays indicating that obscuring 
material (n(HI) ~ 10 22 cm" 2 ) may be surrounding the nuclei. 

We detect X-ray emission coincident with the radio hotspots or lobes in 4 of the 5 targets. This extended emis- 
sion can be explained by the Inverse-Compton (IC) scattering of photons that make up the cosmic microwave 
background (CMB). The magnetic field strengths of ~ 100 — 200 fiG that we derive agree with the equipartition 
magnetic field strengths. The relative ease with which the lobe X-ray emission is detected is a consequence of 
the (1 + z) 4, increase in the energy density of the CMB. For one of the lobes, the X-ray emission could also be 
produced by a reservoir of hot, shocked gas. An HST image of the region around this radio component shows 
bright optical emission reminiscent of a bow-shock. 

By co-adding the 5 fields we created a deep, 100 ksec exposure to search for diffuse X-ray emission from thermal 
intra-cluster gas. We detect no diffuse emission and derive upper limits of ~ 10 44 erg s _1 , thereby ruling out a 
virialized structure of cluster-size scale at z ~ 2. 

The average number of soft X-ray sources in the field surrounding the radio sources is consistent with the number 
density of AGN in the Chandra Deep Fields, with only one of the fields showing a marginally statistically significant 
factor 2 excess of sources with fo.s-2keV > 3 x 10 -15 erg s _1 cm -2 . Analysis of the angular distribution of the 
field sources shows no evi dence for large-scale str ucture associated with the radio galaxies, as was observed in the 
case of PKS 1138-262 bv lPentericci et al] J20021) . 
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1. Introduction 

Radio galaxies can be used to trace the formation and 
evolution of the most massive galaxies known at high red- 
shift. They usually have continuum mor phologies sugges- 
tive of the merging of L* systems (e.g. IPentericci et all 
ll999H200llh and a re surrounded by larg e reservoirs of line 
emitting gas (e.g. Ivan Oiik et al.lll997(l that are compa- 
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rable in size to cD galaxy envelopes. High redshift radio 
galaxy (HzRG) fields targeted by broad and narrow band 
imaging and spectroscopy have been found to locate 'pro- 
toclusters'. They have large excesses of Lya and/or Ha 
emit ters, Lyman break galaxies and extremely red objects 
f fi .g. IPentericci et al.ll2000aUKurk et al .lEoOfl lKurkll200.^ 
IVenemans etaLll2002t Ikilev et al.ll2004h . Thus, HzRGs 
act as beacons to th e progenitors of present-day galaxy 
clusters (see also e.g. IWindhorst et al.lll998l llvison et alJ 
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2000| iBrand et all 120031: ISmail et alJl2003t IStevens et all 
2003). Finding these distant (i.e. z > 2) protoclusters is 



important for constraining models of structure formation 
and cosmology. 

Radio continuum polarimetric observations of HzRGs 
have revealed large rotation measures (RM) and large RM 
gradients. At z > 2 the typical RM observed is several 
hundreds rad m~ 2 , but extreme values of > 1000 rad 
m~ 2 are not uncommon llCarilli et alll997llAthreva et alJ 
ll99SHPentericci et all2000bl) . It is believed that RM arises 
from 'Faraday screens' in our galaxy (i.e. the interstel- 
lar medium, ISM), in intervening galaxies and clusters, or 
near the source itself. At low redshifts, most radio galax- 
ies have relatively small RM values that can be attributed 
to the local ISM CSi mard-Normandin fc Kronberdri980t 
lLeahvlll987L and references therein). However, RM values 
observed towards HzRGs are significantly larger, and are 
believed to be due to Faraday screening of extra-galactic 
origin somewhere along the line of sight. lAthreva et alJ 
l)l998l) show that the chance superpositions of distant ra- 
dio sources and intervening systems such as foreground 
Abell clusters and damped Lya absorbers are too small 
to explain the fraction of high RMs observed. Likewise, 
the expected RMs from such intervening systems are too 
small to account for the high values seen towards dis- 
tant radio sources. Hence, the large RM values may occur 
in a medium close to the source, although it is unclear 
whether they are caused by extended cluster-sized media 
or smaller gas dis tributions on scal es comparable to the 
radio sources fe.g. lCarilli et al.l2002l) . Some radio galaxies 
at low redshift have extremely large values of RM that are 
caused by Faraday rotation in so-called 'cooling flow' c lus- 
ter atmospheres (e.g. 3 C 295; IPerlev fc Ta"vlorlll99lL sec 
ICarilli fc Tavlorl l|2002f) for a review on cluster magnetic 
fields) . There is evidence that the strength of cooling flows 
is correlated with RM, indicating that RM is an indica- 
tor of the fields i n the intracluster medium rather than 
the radio cocoon llTavlor et al.lll994t lEilek fc Owenll2002t 
lEnsslin et al"ll2003l but see lRudnick fc Blundelll l|2003^ . 

X-ray observations are also a powerful tool to find 
deep, gravitational potential wells in the early Universe. 
Diffuse X-ray emission caused by thermal bremsstrahlung 
from a hot ICM, is evidence for a bound system in which 
the gas is in dynamical equilibrium with the galaxies 
and the cold dark matter. Deep X-ray surveys have un- 
covered significa nt numbers of rich galaxy clusters out 
to z ~ 1.2 fe.g. IRosati et al.lll999t IStanford et"ai1l200ll 
iRosati et"al 12004) . Remarkably, the diffuse ICM in these 
distant objects traces the galaxy distribution, and their 
surface brightness and temperature profiles are similar to 
those of lower redshift clusters, indicating that clusters 
formed very early in the history of the Universe. The study 
of thermal emission from HzRGs is important because of 
their suggested linkage to cluster formation. However, if an 
ICM is present in such structures at z > 1.5, the extreme 
cosmological surface brightness dimming makes their de- 
tection very difficult. 



So far only a handf ul of distant radio sou rces has been 
studied with Chandra. iFabian et al.1 |2003) performed a 
very deep (200 ks) study of the radio galaxy 3C 294 at z = 
1.79. They observed a 100 kpc region of diffuse emission 
bounded by sharp edges. Although some thermal compo- 
nent from the ICM could not be ruled out, most of the 
emission was ascribed to Inverse Compton (IC) scattering 
of the cosmic microwave background (CMB) by an older 
population of electrons tracin g out an hourglass- shaped 
region around the radio source . iBelsole et ail l|2004l) found 
weak cluster luminosities a nd IC scatter i ng am ong a sam- 
ple of three bright HzRGs. ICarilli et alJ l|2002ft found dif- 
fuse X-ray emission around the radio galaxy PKS 1138- 
262 at z = 2.16, which is believed to be t he forming, mas- 
sive galaxy at the center of a protoclus ter ( Pentericci et alJ 
l2000al iKurk et aljEfiorl iKurkl 120031) . However, the ex- 
tended emission is seen only along the radio axis, and is 
therefore believed to be associated with shoc ked material 
inside the radio source. The upper limit that lCarilli et all 
(2002) derive for the X-ray luminosity of the ICM (~ 
40% of the l umino sity of the Cygnus A cluster, see e.g. 
ISmith et al.1 l|2002j) ') is not unexpectedly low because there 
has not been enough time since the big bang for these 
struc tures at z ~ 2 to f orm a sufficiently deep potential 
well. Scharf et alJ l(2003h detected extended X-ray emis- 
sion around 4C 41.17 at z = 3.8 roughly following the ra- 
dio morphology. They conclude that the X-ray emission 
arises from 'Inverse-Compton scattering of far-infrared 
photons from a relativistic electron population probably 
associated with past and current activity from the cen- 
tral object', in addition to a lesser contribution from the 
up-scattering of CMB photons. 

It is clear that the nature of the X-ray emission mecha- 
nisms operating in radio galaxies is diverse. Better sample 
statistics are required to enlarge our current understand- 
ing of e.g. the dominant X-ray emission mechanisms, the 
role of magnetic fields on galaxy and cluster-size scales, the 
origin of the hot gas observed in X-ray luminous clusters, 
and large-scale structure associated with HzRGs. Here we 
present observations with the Chandra X-ray observatory 
of five additional radio galaxies at z ~ 2. The sources were 
selected from a compendium of z > 2 HzRGs, focusing on 
the lower end of the redshift distribution in order to min- 
imize surface brightness dimming effects, while still being 
at high enough redshift to allow for cosmological evolution 
effects in source properties. Sources were further selected 
on the presence of one or more of the various character- 
istics: distorted radio morphology, large rotation measure 
(> 1000 rad m~ 2 ), evidence for interaction between Lya 
and radio structure, and extended Lya absorption fea- 
tures. 

Throughout this paper we assume Ho = 70 km s -1 , 
f^M = 0.3, and fl\ = 0.7. The resulting scale factors range 
from 8.0 to 8.4 kpc arcsec -1 from z — 2.6 to z = 2. We 
will use power-law spectra defined as S u ex v~ a , where a 
is the spectral index. 
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Table 1. Log of o bservations. The last colu mn lists the average galactic HI column density in the direction of each 
source as given bv lDickev fe Lockmanl l)l990|) . 



Source 


R.A. (J2000) 


Decl. (J2000) 


z 


Dates of observation 


Exposure time (s) 


n(HI) (xlO™ cm" 2 ) 


0156-252 


0l"58 m 33.5 s 


-24°59'32.26" 


2.09 


2002-11-09 


11477 


1.34 










2002-12-20 


8401 




0406-244 


04 h 08 m 51.5 s 


-24° 18' 16.8" 


2.44 


2002-12-07 


20179 


3.30 


0828+193 


08 h 30 m 53.4 s 


+19° 13' 15.7" 


2.57 


2002-11-05 


19212 


3.72 


2036-254 


20^39™ 24.5 s 


-25° 14' 30.4" 


2.00 


2002-09-01 


19660 


4.92 


2048-272 


20 h 51 m 03.6 a 


-27°03'02.1" 


2.06 


2002-11-11 


19555 


6.29 



2. Observations 

2.1. Observations and data reduction 

Each of the 5 sources was observed for 20 ks with the back- 
illuminated ACIS-S3 chip on Chandra, using the standard 
3.2 s readout timed exposure mode and the faint teleme- 
try format. Sources were positioned at approximately 20" 
away from the aimpoint. The data were processed and 
analysed in June 2003 using the Chandra data analysis 
package CIAO 2.3 1 together with CALDB 2.22. We ap- 
plied the 'ACISABS' script that corrects the Chandra aux- 
iliary response files for a continuous degredation in the 
ACIS quantum efficiency We did not apply the 'tgain' 
correction that has become standard practice only in the 
most recent release of CIAO to correct for the drift in the 
effective detector gain due to an increase in the charge 
transfer inefficiency over time. However, for ACIS-S3 the 
resulting effect on the measured photon energies was of 
the order of only ~ 0.3% near the end of 2002 when the 
data were taken. 

A log of the observations is given in Table 1. 0156- 
252 was observed on 2 separate dates, because the initial 
observation was interrupted after 11.5 ks by a high radia- 
tion shutdown of the system due to an excess of solar wind 
particles. To ease the analysis of these observations, the 
two exposures were merged into a single event file by 're- 
projecting' the events. A lightcurve analysis of the other 
sources indicates no further periods of significant fluctu- 
ations in the count rates of source and large, source-free 
background regions. 

Radio images at 5 and 8 GHz where ob tained by 
ICarilli et all l|l997t) and lPentericci et alJ l|2000ah using the 
VLA in A configuration. The noise is 50 //Jy/beam at 5 
Ghz and 25 /iJy/beam at 8 GHz. The resolution of the ob- 
servations is 0.43" for the 5 GHz maps and 0.23" for the 8 
GHz maps. Analysis of the radio data was performed using 
the Astronomical Image Processing System (AIPS) and 
the Multichannel Image Reconstruction, Image Analysis 
and Display (MIRIAD) software. 

We will also make use of observations obtained with 
the Wide-Field Planetary Camera 2 (WFPC2) on HST. 
0156-252 and 2048-272 were observed as part of Program 
8183 (PI: Miley), and an image of 0406-244 was retrieved 



http:/ /cxc.harvard.edu/ciao/ 



from the HST archive (program 8338, PI: Lehnert). 0156- 
252 was observed for 4800 s and 2048-272 for 7200 s, both 
through filter F555W. This filter has a central wavelength 
of 5443 A. 0406-244 was observed for 2000 s through 
filter F675W which has a central wavelength of 6718 
A. Measurements from these images were made using 
the conversion from counts to flux density as given in 
the PHOTFLAM header keyword of each image. To get 
to the flux density S u (erg s _1 cm -2 Hz -1 ) we calcu- 
late S„ = (X 2 /c) x ( counts /t exp ) x PHOTFLAM (see 
IPentericci et al.lfl999l for a description of procedures). 

2.2. Image registration 

To recover the inherent resolution of the Chandra mir- 
ror/detector system, we removed the pixel randomization 
added in the pipeline processing and regridded our maps 
to 1/10 native ACIS pixel size. After applying a suitable 
smoothing function, we were able to find the position of 
the X-ray core to an accuracy of better than 0.1" pro- 
vided the core had sufficient signal to noise (S/N). We 
then shifted the images to align the X-ray core with the 
radio core. The typical shift needed was one half of an 
original Chandra pixel in right ascension and/or declina- 
tion, consistent with the known astrometric accuracy. In 
the case of 2048-272 no shift was applied, because the nu- 
cleus was not detected in either the radio or the X-rays. 

2.3. Analysis 

The Chandra images are shown in Fig.^ whe re yellow con- 
tours indic ate the 4.7 GHz radio m aps from lCarilli et all 
l|l997|) and lPentericci et all l|2000bj) . We measured the X- 
ray counts from the distinct radio components (e.g. cores 
and lobes) using boxed or circular extraction regions. The 
counts are background-subtracted by using a large, point 
source free, annular region centered on the radio source. 
We restrict our measurements to counts in the energy 
range of 0.2-6 keV for which the contribution of the back- 
ground is minimal. We calculate the total net (i.e. back- 
ground subtracted) counts N = S — B, where S and B are 
the counts in the source region As and the background re- 
gion Ab, respectively. We consider a radio component de- 
tected in the X-rays when N is greater than its formal er- 
ror, a N = [(a s ) 2 + (a B ) 2 } 1/2 , where a s = 1 + (S + 0.75) 1 / 2 



4 



R. A. Overzier et al.: Chandra observations of z ~ 2 radio galaxies 



Table 2. Rotation measures, KMi nt — RM f, s x (1 + z) 2 , 
for the Northern (N.) and Southern (S.) lobes of each 
source taken from the literature. 



Source 


RMi„i 


(rad m 2 ) 


Ref. 




N. lobe 


S. lobe 




0156-252 


1528 ± 296 




1 


0406-244 


880 ± 19 


-705 ± 30 


1 


0828+193 






2 


2036-254 


3 ±227 


-3321 ± 429 


1 


2048-272 


590 




3 



1 Athreya et al. (1998) 

2 Carilli et al. (1997) 

3 Pentericci et al. (2000) 



andas = [1 + (B + Q.75)^ 2 }(A S /A B ) from ICenrelsj l|l 9861) . 
For undetected components we calculate 2a upper limits 
in the 0.2-6 keV band using the above expression for er/v- 
The upper limit on flux is estimated by assuming a power- 
law spectrum with a spectral index of 0.8 using the online 
tool PIMMS 2 . 

The procedure to convert counts into fluxes is as fol- 
lows: we create exposure maps for the soft (0.5-2 keV) and 
hard (2-6 keV) band at nominal energies, and divide the 
images in each band by the corresponding exposure map 
to obtain fluxed images. When extracting the flux for a 
given region, we calculate weight factors defined by the 
ratio of the average energy of counts in that region to the 
nominal energy of each exposure map. The total flux in a 
region is then calculated by summing the weighted fluxes 
in the two bands. We remark that our method is slightly 
different compared to the usual procedure of calculating- 
fluxes from exposure maps using a conversion factor that 
is based on the spectral response within the given band 
as well as an assumed spectral energy distribution. Our 
method is preferred because it allows us to measure fluxes, 
without having to assume a specific spectral shape. 

The measurements are given in Table [3] Where indi- 
cated, fluxes are corrected for galactic absorption using 
PIMMS, assuming a power-law spectrum with a spectral 
index of 0.8 and taking th e HI column density in the di- 
rection of each source from lDickev fc Lockmanl l|l990t in- 
dicated in Tabled All energy ranges are in the observed 
frame, unless stated otherwise. 

3. Results 

3.1. Source description 

We will briefly summarize some of the main char- 
acteristics of each source. Details on radio observa- 
tions, ground-based imaging and spectroscopy, and 
Hubble Space Telescope (HS T) ob ser vations can 
be fo und in iMcCarthv et alJ fl996h. ICarilli et al. 
l|l997f) . IPentericci et all Jl999jK IPentericci et all 

2 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



ll200ll) . | Iwamuro et al.l |2003) and references therein. 
Table [5] summarizes the intrinsic rotation measures, 
RM int = RM obs x (1 + zf. 

• 0156-252 at z = 2.09 has several characteristics rem- 
iniscent of the we l l-studied radio galaxy PKS 1138-262 
(|Carilli et al.ll2002l:|Pentericci et alJl2000a(l . Near-infrared 
continuum shows a host galaxy extended over ~ 2" par- 
allel to the radio axis, and narrow-band Lya shows emis- 
sion line gas extended over the entire region encompassed 
by the northern and southern radio hot spots. The Lya 
brightness is at a minimum at the position of the host 
galaxy and peaks at the location where there is a sharp 
bend in the (northern) radio jet, probably due to shock- 
induced ionization of the gas where the radio jet is de- 
flected by denser material. The northern hot spot has an 
intrinsic RM of ~ 1000 rad m" 2 . 

• 0406-244 at z — 2.44 consists of two main opti- 
cal /infr ared continuum fe atures aligned with the radio 
axis fsee lRush et alJ l|l997fl for a detailed study). The con- 
tinuum components lie embedded in a large Lya halo with 
a spatial extent of ~ 8" along the radio axis. Both radio 
lobes have an intrinsic RM of ~ 1000 rad m~ 2 with oppo- 
site signs. 

• 0828+193 at z = 2.57 is the largest radio source in 
our sample (~ 100 kpc). It has a jet-like feature north 
of the radio core. The northern lobe contains several hot 
spots with a peculiar 90 degree bend. The sou thern radio 
hotspot is unresolved. IPentericci et al ] lll999ll have iden- 
tified the radio core with the brightest optical component 
of several clumps in a ~ 2" region seen by HST. 

No RM is given for this source in Table El since 
the fractional polarization at 8.2 GHz for both lobes 
was below the 4a sensitivity limit of the observations 
(ICarilli et al.lll997t) . 

• 2036-254 at z = 2.00 is the only source in our sample 
for which no optical/near- infrared data are available. The 
radio core is most likely the faint, compact component 
close to the southern lobe. There is a pair of hot spots 
parallel to the radio axis in the north. The southern lobe 
has a RM mt of ~ 3500 rad m -2 . 

• 2048-272 at z = 2.06 is among t he ~ 30% of sources 
in the sample of ICarilli et alJ l|l997|) for which the core is 
undetected. However, a bright near-infrared object seen 
in between the lobes is presumed to be the galaxy hosting 
the radio source. The northern hot spot has a RM ijlt of 
~ 590 rad m~ 2 . 



3.2. X-ray Cores 

We have detected X-ray cores in 4 of the 5 targets. No 
core was detected in the case of 2048-272, and we note 
that the core was also not detected in the radio maps 
indicating that nuclear activity may have ceased. The 
core X-ray luminosities are in the range £x[o.5-6fceV] = 
1.3 x 10 44 - 2.3 x 10 45 erg s -1 . 0156-252 has > 4x the 
luminosity of the other cores. We place an upper limit 
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Fig. 1. Chandra images of 0156-252, 0406-244, 0828+193, 2036-254 an d 2048-272, showing the 0.2-6 keV X-ra y 
images in colorscales, and the 4.7 GHz VLA radio contours in yellow ffrom lCarilli et aDll997trPentericci et alJl2000aj) . 
The 20" x 20" Chandra images have been smoothed using a Gaussian kernel of width 2" (FWHM). Radio contours 
are 1.0 x (0.2, 0.4, 0.8, 1.6, ...) mJy beam" 1 . 
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Table 3. X-ray measurements. 



Source 


e e 


Counts" 


Soft" 


Hard" 


J0.5-2keV 


/2-6/cev 




l°g(-^0.5-6fcev) 




(") 


(0.2-6 keV) 


(0.5-2 keV) 


(2-6 keV) 


erg s _1 cm -2 


erg s _1 cm -2 




erg s _1 


0156-252 


















Total 


4.2 


275 ± 18 


200 ± 14.2 


69 ±8.3 


(3.1 ±0.2) x 10~ 14 


(4.9 ±0.6) x 10" 


-14 


45.42 


Core 


1.1 


249 ± 17 


181 ± 13.5 


62 ±7.9 


(2.8 ±0.2) x 10~ 14 


(4.4 ±0.6) x 10" 


-14 


45.37 


NE 


1.6 


7±4 


4.9 ±2.2 


1.0 ± 1.0 


(7.4 ±3.4) x 10~ ie 


(3.9 ±4.7) x 10" 


-16 


43.56 


SW 


1.2 


4±3 


1.9 ±1.4 


1.9 ± 1.4 


(4.6 ±3.4) x 10" 16 


(1.6 ± 1.2) x 10" 


-15 


43.82 


0406-244 


















Total 


4.8 


21 ±6 


18 ±4.4 


0.5 ± 1.4 


(2.9 ±0.7) x 10~ 15 


(3.6 ±9.3) x 10" 


-16 


44.19 


Core 


0.9 


10 ±4 


8.0 ±2.8 


2.0 ± 1.4 


(1.4 ±0.5) x 10~ 15 


(1.3 ±0.9) x 10" 


-15 


44.11 


SE 


1.8 


9±4 


6.9 ±2.6 


1.0 ± 1.0 


(1.0 ±0.4) x 10" 15 


(5.2 ±6.5) x 10" 


- 16 


43.88 


NW 


2.0 


d 






< 2.4 x 10" 15 




< 44.05 


0828+193 


















Total 


7.0 


32 ±7 


17 ±4.4 


13.3 ±4.0 


(3.2 ±0.8) x 10~ 15 


(1.2 ±0.4) x 10" 


- 14 


44.91 


Core 


1.0 


23 ±6 


12 ±3.5 


9.9 ±3.2 


(2.1 ±0.6) x 10 -15 


(8.7 ±2.8) x 10" 


-15 


44.78 


NE 


1.4 


d 






< 1.9 


x 10~ 15 




< 44.00 


SW 


1.0 


d 






< 1.6 x 10" 15 




< 43.93 


2036-254 


















Total 


4.1 


35 ±7 


14 ± 3.9 


18.5 ±4.4 


(2.6 ±0.7) x 10~ 15 


(1.2 ±0.3) x 10" 


-14 


44.65 


Core 


0.6 


22 ±6 


7.0 ±2.6 


15.0 ±3.9 


(1.4 ±0.5) x 10~ 15 


(9.6 ±2.5) x 10" 


-15 


44.51 


NE 


1.3 








< 2.0 x 10~ 15 




< 43.76 


SW 


1.4 


6±4 


5.9 ±2.5 


0.0 ±0.0 


(1.0 ±0.4) x 10 -15 






43.47 


2048-272 


















Total 


4.3 


8±4 


4.1 ±2.2 


4.0 ±2.2 


(6.0 ±3.3) x 10~ 16 


(2.9 ± 1.6) x 10" 


■15 


44.05 


Core 


0.8 


d 






< 2.1 x 10" 15 




< 43.82 


NE 


1.7 


6±4 


2.8 ± 1.7 


1.8 ± 1.4 


(4.5 ±2.8) x 10~ 16 


(1.2 ± 1.0) x 10" 


-15 


43.73 


SW 


1.5 


d 






< 1.7 x 10 -15 




< 43.72 



a Errors are calculated using Gehrels (1986) in the low-count regime. 
b Errors are V counts. 

c Observed fluxes (i.e. not corrected for galactic absorption). 

d Undetected. We give a 2a upper limit for the total flux at 0.2-6 keV. 

e Circular extraction radius. 



on the core luminosity of 2 048- 272 of 7 x 10 43 erg s" 1 . 
lHardcastle fc Worralll (|l999j) and lBrinkmann et all (|2000() 
presented a correlation between the X-ray and radio lu- 
minosity of radio-loud quasars from the cross-correlation 
of ROSAT and radio sources. It has been suggested that 
this correlation implies a physical relationship between 
X-ray and radio cores, which can be explained by mod- 
els in which e.g. the X-ray emission originat es from the 
base of the jet feiardcastle fc Worralll Il999l). In Fig. 
we show objects from the iBrinkmann et alJl)20f)0|) sam- 
ple classified as narrow-line AGN, as well as the best- 
fit correlation for their subsample of radio-loud quasars, 
\ogLx — 11.22 ± 0.483 xlogL r . We have also indicated 
several other z > 2 radio galaxies (1138-262 at z — 2.16, 
B0902±343 at z = 3.395, 4C41.17 at z = 3. 798) and the 
local source Cygnus A for comparison (se e ICarilli et alJ 
l2002tlScharf et alJ2003HFabian et all2002l and references 
therein). 

As shown in Fig. [21 three out of four sources in our 
sample agree with the L r /Lx relation if we use the ob- 
served 2-6 keV (rest-frame 6-18 keV) flux with a power- 
law spectrum (a = 0.8) to calculate the luminosity at rest- 



frame 2 keV (observed 0.7 keV). This result is virtually in- 
dependent of the amount of intrinsic absorption, since it 
is of negligible effect at rest-frame 6-18 keV. Calculating 
the 2 keV luminosity from the 0.5-2 keV observed flux 
(rest-frame 1.5-6 keV) yields values that are significantly 
lower, but for 0156-252, 0828±193 and 2036-254 this can 
be explained by an intrinsic absorption column density of 
rc(HI) ~ 10 23 cm" 2 . Interestingly, the core of 0406-244 
seems highly underluminous in Fig. a discrepancy that 
cannot be explained away by invoking a large amount of 
absorption. However, this may be a S/N issue: we have de- 
tected only 2 photons in the hard 2-6 keV band, compared 
to 62, 10 and 15 for the other 3 sources. If we instead cal- 
culate the 2 keV luminosity from the 0.5-2 keV observed 
flux (8 photons detected) with n(HI) = 10 23 cm -2 we 
find \og{L2kev) ~ 26.4 compared to ~ 25.9 for the ear- 
lier method, and note that its deviation from the L r /Lx 
relation is now comparable to that of 4C41.17. 

0156-252 is the only source in our sample for which 
the core is detected with sufficient S/N to carry out a 
(crude) spectral analysis. The X-ray spectrum is shown in 
Fig. |21 Assuming a power-law spectrum with absorption 
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Fig. 2. Monochromatic X-ray luminosity versus radio luminosity. The four detected cores from our sample are indi- 
cated by filled stars. Their X-ray luminosities were corrected for galactic absorption only. Narrow line AGN from the 
FIRST/ROSAT sample of iBrinkmann et alJ feOOd) are indicated by circles. The dashed line indicates the correlation 
for radio-loud quasars of lBrinkmtinn*'et alJ 1 2000^ Open stars indicate several low (Cygnus A at z — 0.06) and high 
reds hift (1138-262 at z = 2.16, B0902+343 at z = 3 395, 4 C41.17 at z = 3.798) radio galaxies from the literature 
fsee lCarilli et al.ll2002t IScharf et al1l2003t iFabian et alJl2002l and references therein). The X-ray luminosities of these 
sources were all corrected for intrinsic absorption. 



3.3. X-ray emission from lobes/hotspots 

We obtained positive detections of X-ray components co- 
incident with the northeastern (NE) and southwestern 
(SW) radio lobes of 0156-252, the southestern (SE) lobe 
of 0406-244, the SW lobe of 2036-254, and the NE lobe 
of 2048-272. Below we will evaluate the likely mechanisms 



> 



at the redshift of the source, we find ax = 0.8 ± 0.2 and 
n(HI) = (1.6 ± 0.7) x 10 22 cm" 2 . The fit has a reduced \ 2 
of 0.34, and is indicated in Fig. EI The error bars for the fit 
parameters correspond to a change of 1 in the reduced y 2 ■ 
Repeating the fit using using the Cash statistic valid for ^ 
the low counts regime we find n(HI) = (1.5 ± 0.8) x 10 22 » 
cm~ 2 with a unchanged. The galactic HI column density J2 
towards this source is 1.34 x 10 20 cm -2 (see Table 1). For 3 
completeness, we note that the spectrum can be fitted by 
a Mewe-Kaastra-Liedahl (Mekal) thermal plasma model 
with kT = 23 ± 8 keV (reduced \ 2 of 0.50). 

The signal-to-noise of the 3 other cores is insufficient 
to make spectra. However, if we assume that they are a 
generic class we might estimate the amount of intrinsic 
absorption by stacking the data. The hardness ratio HR = 
(H — S)/ (H + S) (where H and S are the number counts 
observed in the hard and soft bands) is 0.0±0.2, consistent 
with a Type II AGN at z ~ 2 with an obscured (n(HI) « 
10 23 cm " 2 ) power-law spe ctrum with a photon index of 
r = 1.8 (|Tozzi et al.ll200lh . 



0.01 - 



o 0.001 - 




2 3 4 

Energy (keV) 

Fig. 3. Spectrum of the core of 0156-252, fitted with a 
power-law spectrum modified by absorption at the source. 



for these detections. Given the generally low signal to noise 
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of all features we detected, we will focus on the X-ray emis- 
sion that is identified with the radio lobes 3 . 

• Synchrotron radiation 

In Fig.0]we plot spectra for the lobes that have a detection 
in the X-ray. The spectra were constructed by (i) extrap- 
olating the radio flux densities to X-ray frequencies using 
the radio spectral index af; 2 -, and (ii) using the 4.7 GHz to 
X-ray spectral index, af 7 . The results of our analysis are 
indicated in Tabled] In the first method (indicated by dot- 
ted lines in Fig. [3}, we can directly compare the observed 
X-ray flux density to the predicted synchrotron flux den- 
sity based on the radio spectrum. Except for the NE lobe 
of 0156-252, the predicted flux densities are several orders 
of magnitude lower than that observed. This seems to rule 
out a synchrotron origin for the X-ray emission from these 
lobes. However, S v is a strong function of the spectral in- 
dex, since we are extrapolating over ~ 10 s decades of fre- 
quency. Only a small offset in af% can lead to a significant 
under- or overprediction of the X-ray flux density. The ra- 
dio spectral indices used were measured over a relatively 
narrow frequency range, and may be subject to errors. 
We estimate that the typical error in af;f is ~ 0.2, based 
on the comparison with a spectral index a\ \ determined 
by combining our data with data (not shown here) from 
the publicly available NRAO VLA Sky Survey (NVSS) 4 , 
albeit at much lower resolution. Therefore, we attempt 
the second method (indicated by solid lines in Fig.0J. We 
slightly modify the above procedure and now estimate the 
radio flux density expected at 8.2 GHz, assuming a priori 
that the spectrum is synchrotron with a radio to X-ray 
spectral index a* 7 . The ratio of predicted to observed 8.2 
GHz flux densities is 1.0 for the NE lobe of 0156-252 and 
ranges from 1.2 to 1.5 for the other lobes. This indicates 
that the observed X-ray flux densities are consistent with 
synchrotron radiation, provided that the error in the in- 
tegrated flux densities at 8.2 GHz are > 20%. Typically, 
errors in the integrated radio flux densities are of the order 
of ~ 5%. Furthermore, the extrapolation from the radio to 
the X-ray using the radio spectral index wnrier-predicted 
the observed X-ray emission in all cases. If there were to 
be a large, random error in the integrated radio flux den- 
sities, we might have expected to see some of the spectra 
ouer-predict the X-ray flux as well. The fact that this is 
not observed may indicate that the typical error in the 
radio flux densities is indeed only a few percent. 

To further test our synchrotron models, we have ex- 
tracted flux densities in the optical (~ 6000A) using the 
HST/WFPC2 observations of 0156-252, 0406-244, and 
2048-272 (see Sect. EJ. In Fig. |we have indicated la 
upper limits on the optical emission of the SW lobe of 
0156-252, the SE lobe of 0406-244 and the NE lobe of 
2048-272. Although the flux densities in the optical are 



3 Our radio data are inadequate to confidently identify radio 
features as lobes or hotspots. We adopt the term lobes to mean 
any radio emission not associated with the core. 

4 http://www.cv.nrao.edu/nvss/ 



consistent with the synchrotron models, they are unfortu- 
nately not deep enough to rule them out. 

Interestingly, the optical flux density of the NE lobe of 
0156-252 is about 1 order of magnitude higher than pre- 
dicted. We suspect that the optical emission in this lobe 
is most likely not continuum emission, but line emission. 
The NE lobe coincides with the peak of the Lya intensity, 
and several strong emission lines common to HzRG spec- 
tra fall within the filter. Moreover, the HST/WFPC2 im- 
age (Fig. EJ) shows a bright component roughly following 
the morphology of the radio lobe, reminiscent of a shell 
of shocked gas. We will treat this feature in a separate 
discussion below. 

• Synchrotron self-Compton (SSC) emission 
SSC emission arises due to Inverse- Compton scattering of 
local synchrotron photons. This process dominates over 
other IC mechanisms when the local synchrotron pho- 
ton energy density is higher than that of the external 
(e.g. CMB) photon field, and has been found to explain 
the Chandra hot s pots of radio sou rces such as 3C 29 5 
l)Harris et al.ll200(t and Cygnus A (jWilson et a.l J l200ftl . 
We calculate the energy density of the local synchrotron 
photon field, u s , by integrating the radio flux density over 
1 decade of frequency (1-10 GHz) and assuming a cylin- 
drical geometry for the radio lobes. The results are listed 
in Table El 

A detailed calcul ation of SSC (see e.g. 
iBand fc Grindlavl lll985h for the theoretical frame- 
work) would require knowledge of the full radio spectrum 
and the geometry of lobes and hot spots. However, we can 
roughly estimate the SSC flux density in the X-ray from 
the ratio of the energy losses in the IC and synchrotron 
channels H=u s /ub ~ Lic/Ls, where ug is the energy 
density of the magnetic field B, and Ljn and L,g are the 
IC and synchrotron luminosities fsee lHarris et alJ [2000: 
Donahue et alJl2003|) . We calculate ug from the magnetic 
field strength B eq that will make the energy densities 
of fields and par ticles a pprox imately equal. We use the 
formula given bv iMilevI (|1980): 

B eq = 5.69 x 10~ 5 x 

2 

" (1 + fc) (1 + zf +a -S r v^ v\~ aT - tj~° r l 7 
V 9 2 s sin^ 4> \ ~ ot r 

where B eq is in Gauss, k is the ratio of energy in heavy par- 
ticles to that in electrons, ij the filling factor, 9 2 (arcsec 2 ) 
the area taken up by the radio lobe (assumed to be a 
cylinder viewed broadside), s the pathlength through the 
lobe (kpc), and v\ and (GHz) are the lower and up- 
per cut-off frequencies of the synchrotron spectrum. We 
take = 7r/2,fc = O and filling factor unity so as to 
obtain a real minimum, and note that the results are 
only weakly dependent on these parameters given the 2/7 
power. R, B eq and Sx,ssc are given in Table [5] where 
B eq ~ 100 — 200 /zG. In all cases, the ratio of predicted 
X-ray SSC flux to observed X-ray flux is less than 0.4. 
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Fig. 4. Synchrotron model predictions. We indicate the observed radio, optical and X-ray flux densities for all com- 
ponents that have a positive X-ray detection. Solid lines represent synchrotron spectra constructed using the 4.7 GHz 
radio to X-ray spectral index, af 7 . Dotted lines indicate spectra constructed using a radio spectral index, a^ y. For each 
component, the data and models were offset for clearer visibility of the results. The multiplicative factors for 0156-252 
NE, 0156-252 SW, 0406-244 SE, 2036-254 SE and 2048-272 NE were 1.1, 1.4, 0.03, 0.001, 0.0003, respectively. 

Table 4. Synchrotron models. 



Component 


°A.7GHz 

(mJy) 


"8.2GHz 

(mJy) 


s.ib 

«4.7 


rtc 

° X , model 

(Jy) 


^X,obs 

(Jy) 


Q4.7 


8.2GHz, model 

(mJy) 


S 9 . 

5000A, model 

(Jy) 


0156-252 NE 


89 


47 


1.15 


5.0 x 10" 11 


8.5 x 10" 11 


1.13 


48 


1.6 x 10"' 


SW 


7 


3.5 


1.25 


6.4 x 10" 13 


1.5 x 10" 1() 


0.96 


4.1 


9.2 x 10" 8 


0406-244 SE 


31 


11 


1.86 


4.8 x 10" 1 ' 


1.2 x 10" 1 " 


1.05 


17 


1.4 x 10"' 


NW 


72 


32 


1.46 


1.5 x 10" 13 


< 1.8 x 10" 10 


> 1.07 


< 40 


< 2.4 x 10" 7 


0828+193 NE 


7 


2.0 


2.25 


1.0 x 10" al 


< 1.4 x 10" 1U 


> 0.96 


< 4.1 


< 8.7 x 10 -8 


SW 


9 


3.6 


1.65 


6.0 x 10" 16 


< 1.2 x 10" 10 


> 0.98 


< 5.2 


< 8.6 x 10" 8 


2036-254 NE 


51 


29 


1.01 


3.7 x 10" 1 " 


< 1.5 x 10" iu 


> 1.06 


< 28 


< 1.9 x 10"' 


SW 


37 


15 


1.62 


4.2 x 10" 15 


7.7 x 10" 11 


1.08 


20 


1.1 x 10" 7 


2048-272 NE 


83 


34 


1.60 


1.4 x 10" 14 


1.3 x 10" 1 " 


1.10 


45 


2.0 x 10"' 


SW 


6 


1.1 


3.05 (2.2 used) 


2.1 x 10" 20 


< 1.3 x 10" 10 


> 0.96 


< 3.5 


< 7.7 x 10" 8 



a Integrated radio flux density at 4.7 and 8.2 GHz. The typical error is assumed to be ~ 5%. 

6 Radio spectral index. The typical error is assumed to be ~ 0.2. 

c Predicted flux density at 2 keV for a radio synchrotron model. 

d Observed flux density at 2 keV. 

e X-ray to (4.7 GHz) radio spectral index. 

f Predicted flux density at 8.2 GHz. 

9 Predicted flux density at 5000A. 



• Inverse- Compton (IC) scattering of the CMB 
Another possibility is the up-scattering of CMB photons 
by relativistic electrons in the radio source. This effect is 
expected to become significantly more dominant at higher 



redshifts due to the (1 + z) 4 increase in the energy density 
of the CMB (u CMB ~ 3.4 x 10" 11 erg cm" 3 at z ~ 2). 

The scattering electrons are assumed to belong to the 
same power law distribution of electrons which is respon- 
sible for the radio synchrotron emission. Therefore, the 
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Table 5. Predictions for SSC and IC/CMB emission, and the equipartition and IC/CMB magnetic field strengths. 



Component 




area" 




b 

u~ 


/ c 

Us/UCMB 


H 


re 

J0.5-6keV,SSC 


fo.5— 6keV,SSC ^ 
f().5-6keV,obs 


DO 

B a eq 


D IC/CMB 


0156-252 NE 


r 


= 0'.'35, I = 


= 170 


2.8 


0.7 


0.04 


4.4 


0.4 


129 


104 


SW 


r 


= 0'.'4, I = 


077 


0.3 


0.1 


0.02 


0.2 


0.01 


67 


30 


0406-244 SE 


r 


= 0'.'3, I = 


172 


5.6 


0.9 


0.05 


2.7 


0.2 


169 


179 


NW 


V 


= Of! 3, I = 


172 


6.0 


1.0 


0.05 


5.1 


> 0.2 


170 


> 130 


0828+193 NE 


T 


= 0'.'2, I = 


0'/7 


7.9 


1.2 


0.04 


0.7 


> 0.04 


213 


> 152 


SW 


T 


= 0'.'2, I = 


077 


3.3 


0.5 


0.03 


0.4 


> 0.03 


162 


> 93 


2036-254 NE 


r 


= 073 




2.9 


0.9 


0.04 


2.1 


> 0.1 


139 


> 44 


SW 


r 


= o"3, i = 


170 


2.7 


0.8 


0.04 


2.0 


0.2 


135 


135 


2048-272 NE 


r 


= 075, I = 


171 


3.4 


0.9 


0.04 


6.7 


0.4 


126 


154 


SW 


T 


= 076, I = 


274 


0.3 


0.1 


0.02 


0.3 


> 0.02 


61 


> 117 



a Cylinder of length I and radius r (r along the line of sight), or a sphere of radius r. 

b Synchrotron energy density, its, in units of 10 -11 erg cm -3 . 

c Ratio of the synchrotron to cosmic microwave background energy densities. 

d Ratio of the synchrotron to equipartition magnetic field energy densities, R — u$/ub- 

e Predicted synchrotron self-Compton (SSC) flux in the 0.5-6 keV band in units of 10 -16 erg s _1 cm 

* Ratio of predicted SSC flux to observed X-ray flux. 

9 Equipartition magnetic field strength in fiG. 

h IC/CMB magnetic field strength in fiG. 



X-ray flux may be used to constrain the magnetic field 
strength: since the energy density of the CMB is fixed at 
any given z, less observed X-ray flux implies fewer elec- 
trons, so the magnetic field producing the observed radio 
synch rotron emission must be s tronger. Using the formal- 
ism o f lHarris fc OrindlavHl979|) we calculate the magnetic 
field strength, Bjc'- 

1+ar _ [ (5.05 x 10 4 )^C(a r )G(a r )(l + zf+ a ^S r u^ ' 
IC ~ [ lO^&ti/g" 

where Bic is in Gauss, a r is the radio spectral index, 
S r and Sx are the radio and X-ray flux densities (erg 
s^ 1 cm -2 Hz^ 1 ) at v r and vx (both in Hz), respectively. 
C(a r ), which is well approximated by the value 1.15xl0 31 , 
and G(a r ), which is a function that is slowly y arying with 
a r , can be found in lHarris fc Cxrindlavl il97$ . From the 
radio and X-ray flux densities listed in Table 21 we derive 
field strengths of ~ 30 — 180 /iG for the components that 
are detected in the X-ray. The results are listed in Tabled 
where we also give lower limits on Bjc for the undetected 
lobes. 

We can compare Bjc to B eq , which is an estimate of 
the magnetic field that is solely based on the observed ra- 
dio synchrotron flux (see expression above) . The B eq field 
strengths are typically ~ 100 — 200 /xG, remarkably close 
to the field strengths derived for the IC/CMB mechanism. 
The general agreement that we find between the magnetic 
field strengths using the two independent field estimators 
is consistent with the X-ray flux being produced by the 
IC/CMB process. 

The (observed) 1-10 GHZ synchrotron flux in HzRGs 
is produced by relativistic electrons that have 7 ~ 
10 3.4-3.9 for B _ 10Q ^ syn = 4.2(B/l A1 G)7 2 Hz, 



see lBaechi et al.lll998|) . The number density of relativis- 
tic electrons with energies between 7 and 7 + ^7 can 
be expressed as a power-law, N^dj = N^ s d r y with 
s = 2a + 1. The up-scattering of CMB photons to X- 
ray frequencies is provided by electrons with a Lorentz 
factor 7 ~ 10 3 , since v out m 7 2 z^„, and the frequency 
for which the energy density of the CMB peaks is ~ 
1.6 x 10 11 x (1 + z) Hz. Electrons having such Lorentz 
factors are highly abundant, given the extrapolation of 
N(fy) to 7 ~ 10 3 . Note that in the SSC process described 
above, the X-ray emission must be produced by the up- 
scattering of synchrotron photons off a population of rela- 
tivistic electrons having 7 ~ 10 4 ~ 4 ' 5 . Since Nfa) decreases 
strongly with increasing 7 and a > 1 for all our sources, 
iV(10 4 ' 5 ) will be several orders of magnitude lower than 
7V(10 3 ). Together with the fact that u s /ucmb, is around 
unity (see Table 0) this implies that the IC/CMB process 
will be a far more efficient process for producing X-rays at 
these redshifts. While the contribution of SSC to the ob- 
served X-ray flux might still be in the range of ~ 1 — 40% 
as detailed above, the SSC output is expected to peak 
around 10 15 Hz, or in the UV part of the spectrum. 

• Other mechanisms 

The ambient medium in which a radio source expands is 
expected to be very different for HzRGs compared to radio 
sources at low redshift. Low redshift sources often lie in 
a smooth, virialised (T ~ 10 8 K) atmosphere. For HzRGs 
such a smooth medium does probably not exist because 
the necessary potential well has not yet formed. Instead 
they lie in a multiphase medium, consisting of cold (10 
K), high density clouds embedded in low density regions 
approaching virial temperatures (10 6-7 ). The passage of 
a radio jet through the multiphase can lead to interesting 
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Fig. 5. HST/WFPC2 F555W image of 0156-252 with 4.7 GHz VLA radio contours superimposed. We have detected 
extended UV line emission (most likely CIII] A1909 A) at the position of the NE lobe (box region). An elongated 
feature in the extended emission resembles a bow-shock that follows the radio lobe (arrow). 



phenomena, such as jet-induced star formation in the high 
density clouds and shock heating of lower density regions. 
When intermediate density regions are shocked they may 
cool off, emitting strong emission lines such as Lya. The 
lower density gas with high filling factor will be shock 
heated to temperatures where the gas starts producing 
X-rays. Can the X-ray emission coincident with the NE 
lobe of 0156-252 be produced by shocks in this thermal 
gas? Interestingly, the NE radio lobe of 0156-252 coincides 
with the peak of the Lya emission, and the HST/WFPC2 
image (see Fig. |3J) shows evidence for a shell of shocked 
(emission-line) ga s that follows the be nd of the radio lobe. 

For 1138-262. ICarilli et alJ l)2002|) propose that much 
of the extended X-ray emission comes from ambient gas 
that is shock heated by the expanding radio source. For 
this gas they estimate a density of 0.05 cm -3 and a pres- 
sure of 10 -9 dyn cm -2 . This pressure is comparable to the 
optical line emitting gas and to minimum pressures in the 
radio source, with a total 

gas mass of about 2.5 x 10 12 Mq. Carilli et al. hy- 
pothesize that the high filling factor X-ray emitting gas 
may confine both the radio source and the line emitting 
clouds. The extended X-ray luminosity associated with 
the NE lobe of 0156-252 is a factor 5 or so less than in 
1138-262, implying a factor 5 lower total gas mass (for a 
fixed density) if the X-rays represent shocked gas. 

So far, the discussion has been limited to the X- 
ray emiss ion directly coin cident with the radio lobes and 
hotspots- lBrunettil ifeoOOh describes a mechanism in which 
IC scattering of (mostly infrared) AGN photons by rela- 
tivistic electrons can also produce X-ray bright emission. 
Such a mechanism, if it exists, may contribute significantly 



to the radio/X-ray alignment effect. Under the Brunetti 
mechanism the X-ray luminosity is expected to decrease 
with increasing distance from the source (i.e. the hidden 
quasar) producing the photons. This will result in opposite 
gradients in the X-ray and radio luminosities, and entails 
(in our case) t he presence of a n inner, undetected part of 
the r adio lobe. lBrunetti l)2000|) predicts that the receding 
radio lobe will produce brighter X-ray emission than the 
approaching lobe, because time delay makes it closer to 
the nucleus and because backward scattering is more ef- 
fective than forward scattering. In 0406-244 (and possibly 
in one or several of the other sources as well) there is X-ray 
emission in between the core and the SE lobe. The nature 
of this emission could be similar to the aligned emission 
seen in several other high-redshift quasars and radio galax- 
ies (e .g. lYuan et al.ll2003t iFabian et ai]l2003l IScharf et all 
2003) for which the Brunetti model is among the possible 
scenarios. However, the proposed mechanism relies on a 
(suggested) supply of electrons with 7 < 100. Given the 
highly speculative amplitudes at these low energies, we 
will not pursue this mechanism quantitatively. 

3.4. Thermal emission from hot (cluster) gas 

An important driver for observing powerful HzRGs in the 
X-ray is to search for traces of thermal emission from an 
ICM. To test whether our data show any evidence for the 
presence of extended, diffuse X-ray emission we will at- 
tempt two methods: 

• Smoothed fields 
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Fig. 6. X-ray image of the field around 0156-252 . The image was smoothed by a Gaussian (10" FWHM), after 
replacing the X-ray counts due to the radio source and point sources in the field by a Poissonian background. A large 
region of diffuse emission (indicated by the dashed ellipse) is found just northeast to the radio galaxy (indicated by 
its radio contours). 



To obtain upper limits on the luminosity of extended 
regions in each field, we place 10 circles of 12" diameter 
(corresponding to ~ 100 kpc at z ~ 2) around each radio 
source in regions without visible point sources. We calcu- 
late the average background count in these regions and 
the lcr deviation from the mean. We use the resulting 5er 
countrates to calculate 5er upper limits on the X-ray flux 
of typical 100 kpc-sized areas. We assumed a Raymond- 
Smith thermal spectrum with kT = 1 keV, and used the 
galactic uhi for each source to produce the unabsorbed 
flux at 0.2-6 keV. We derive an upper limit on the flux of 
(4-9) xlO -15 erg cm -2 s _1 , corresponding to luminosities 
of (1.5-4) xlO erg s _1 (for a uniform sphere). 

A 5er deviation from the background countrate over a 
12" region would easily stand out from the X-ray maps 
after smoothing with a large kernel. Therefore, we cre- 
ate smoothed images of diffuse emission. First, we remove 
counts that are associated with the radio galaxy by re- 
placing all counts in a circular region encompassing the 
entire radio structure by the background using the task 
DMFILTH that maintains the Poissonian nature of the 
background. The background was estimated from a large 
annular region around the radio source. We then run the 
WAVDETECT algorithm (Freeman et al. 2002) within 
CIAO to identify the remaining point sources and replace 
them by the background. The images were smoothed by 
a 10" (FWHM) Gaussian. The images were then visually 
inspected to look for regions of diffuse emission. We found 
a single detection of an extended X-ray component in the 
smoothed map of 0156-252 shown in Fig. [I)] In an ovally 
shaped region that can be approximated by an ellipse of 
10" x 17" (83 x 142 kpc if it were to be located at the same 
distance as the radio galaxy), we find 14.4 ± 5 net counts 
(3<r) in the energy range (0.2-6) keV with an average en- 



ergy of 1.2 keV. The centroid of the region is situated 
approximately 12" to the northeast of the core of the ra- 
dio source. The (observed-frame) 0.5-2 keV unabsorbed 
flux is 2.3 x 10~ 15 erg cm~ 2 s _1 . If this emission were 
to come from a sphere of thermal gas at the redshift of 
the radio source, it would correspond to a luminosity of 
7.4 x 10 43 erg s" 1 and a mass of 10 11 Mq (assuming a 
uniform sphere). We could not identify this X-ray region 
with objects in any of our optical and near-infrared imag- 
ing data. Therefore, in this particular case the true nature 
of the emission remains highly obscure, and we cannot rule 
out it being due to noise. However, it illustrates the kind 
of emission that might be detectable in the search for high 
redshift cluster gas, provided a proper identification can 
be made. 

• Stacked fields 

To further search for emission from thermal gas in an an- 
nular region surrounding the radio sources, we stacked the 
five 20 ks exposures to obtain a single 100 ks exposure of a 
typical HzRG field. Each field was cleared of point sources 
before stacking them. In the stacked image, we then re- 
placed the entire region inside a circle encompassing the 
largest of the radio sources by the background to ensure 
that we are not measuring extended X-ray emission asso- 
ciated with the radio structures. The stacked image was 
smoothed by a 10" (FWHM) Gaussian. The stacked im- 
age and the smoothed image are shown in Fig. The 
region used for the background determination is indicated 
by the large, dashed annulus. The small circle in the center 
of the field indicates the size of the largest radio source. 
From the smoothed image, where the scales run from 2cr 
below the mean to 2a above the mean, we see no evidence 
for diffuse, extended emission in the vicinity of the radio 
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Fig. 7. Left: Stacked image of the five HzRG fields, after subtraction of point sources. Indicated are the annuli from 
which a radial surface brightness profile was extracted. The innermost circle encompasses the total extent of the largest 
radio source in the sample. Counts inside this region due to the radio sources have been replaced by the background, 
estimated from the large, dashed annular region. Right: Stacked image smoothed by a 10"(FWHM) Gaussian. Scales 
run from 2a below (white) to 2a above (black) the mean. The background region and the maximum size of the radio 
sources are indicated. The distinct peaks in this summed image are not traced back to significant features in any of 
the individual images. 



source(s). We obtained a radial profile of surface bright- 
ness from the stacked image using the annuli as indicated 
in Fig. [7] The radial profile out to 1' is consistent with 
zero contribution from diffuse, extended X-ray emission, 
as shown in Fig. |H1 Out to a radius of ~ 75" (~ 625 kpc) 
we measure 25 ± 85 counts. Assuming a 1 keV thermal 
spectrum and correcting for galactic absorption, the 3a 
upper limit on the (observed) 0.2-6 keV flux inside this 
radius is 1.2 x 10~ 14 erg cm -2 s _1 , corresponding to a lu- 
minosity of < 4x 10 44 erg s _1 . We derive a central electron 
density of < 7 x 1CP 3 cm~ 3 and a total enclosed gas mass 
of < 4 x 1O 13 M0, assuming an isothermal sphere with a 
cluster /3-model surface brightness profile with [3 — 0.67 
and a core radius of 200 kpc. 

3.5. Serendipitous X-ray sources 

IPentericci et alJ l|2002|) found evidence of an excess of 
Chandra detected X-ray sources in the vicinity of radio 
galaxy PKS 1138-262 at z = 2.16. They found 16 sources 
in the soft band (0.5-2 keV) with a minimum flux of 
1 x 10~ 15 erg cm~ 2 s" 1 , compared to 10.3 (11) calculated 
from the Chandra Deep Field South (North) for a simi- 
larly sized region, and 8 sources with fluxes > 3 x 10 -15 
erg cm~ 2 s _1 compared to 4 (5.1) for the Chandra Deep 
Fields. Although this may not be considered a substan- 
tial overdensity in angular space (an excess of ~ 50% 
compared to a typical cosmic variance of 20 — 30%), it 
is very significant in redshift space given that six of the 
serendipitous X-ray sources could be identified with pre- 
viously discovered Lya or Ha emitting galaxies in the 
protocluster surrounding PKS 1138-262 l)Pentericci et alJ 



2002). Therefore, observing overdensities of X-ray sources 
may provide evidence for the existence of gal axy overden- 
sities associated with radio sources (see also ICappi et alJ 

Soo3). 

Following IPentericci et al 

] |2002), we use 
WAVDETECT to investigate the presence of serendip- 
itous sources in our 5 fields. Because the 2 separate 
observations of 0156-252 were taken at different roll 
angles, the corners of the ACIS-S3 chip only have half of 
the total 20 ks exposure time. We specify a circular region 
covering the maximum area having the full 20 ks depth. 
The circle has a radius of 3.852', thereby covering an area 
of 46. 62' 2 . For consistency we use this region in all the 
fields. We limit the ene rgy range for source de tection to a 
maximum of 2 keV (cf . IPentericci et al.ll2002|) to improve 
the S/N (the background is lower at lower energies). 
For each field, we exclude the radio source, and we only 
consider sources that are detected with S/N > 3, where 
we define S/N as the number of background-subtracted 
counts divided by the root of the total counts. The net 
count rates of the detected sources were converted into 
fluxes assuming a power-law spectrum with Galactic 
absorption along the line of sight, and a photon index of 
r = 2.0 for the soft band. Table lists the number of 
sources in each field that have soft (0.5-2 keV) fluxes of 
> 1.5 and > 3 x 10~ 15 erg cm~ 2 s _1 . Besides the number 
of sources found in each 46. 62' 2 area we give the number 
density per square degree between parentheses. We also 
list the number of sources in each bin averaged over the 
5 fields. For comparison, we have run the same detection 
procedure on the observations of PKS 1138-262 using 
(i) the full 40 ks exposure, and (ii) a 20 ks subsample. 
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Fig. 8. Radial surface brightness profile extracted from the stacked image shown in Fig. 

Table 6. The number of X-ray sources detected in the fields of z ~ 2 radio galaxies, and comparison with the field 
population expected based on the Chandra deep fields. 



Field 


Number of sources" 




/o.5-2fceV > 1-5 


/o.5-2fceV > 3-0 


0156 - 252 


7 (541) 


5 (386) 


0406 - 244 


7 (541) 


7 (541) 


0828 + 193 


5 (386) 


4 (309) 


2036 - 254 


7 (541) 


5 (386) 


2048 - 272 


2 (154) 


1 (77) 


average of fields 


6 (463) 


4 (309) 


1138 - 262 c 


7 (541) 


6 (463) 


1138 - 262 d 


11 (715) 


7 (541) 


CDFs e 


6/7 (470/540) 


3/4 (260/330) 



a Number of X-ray sources sources found in the field (per deg 2 ) 
b Minimum source flux in units of 10 -15 erg s _1 cm" 2 

c Comparison 1: A 20 ks subsample of the full 40 ks o bservation of 1138-262 
d Comparison 2: The full 40 ks exposure of 1138-262 JPentericci et alj|2002l) 

e Comparison 3: The Chandra Deep Field South /North jMushotzkv et alJl20o3 ICiacconi et alJl200ll) 



We have also indicated the n umber of sources found 
in the Chandra D eep Fields l)Mushotzkv et all l2000l 
iGiacconi et alEoOlli. 

The number of bright sources in the field of 0406-244 is 
twice as high compared to the Chandra Deep Fields, and 
similarly high as PKS 1138-262. The number of sources 
in the field of 2048-272 is less than half of the number 
of sources in the CDFs in both flux bins. On average we 
find 6 and 4 sources in the > 1.5 and > 3 x 10~ 15 erg 
cm~ 2 s _1 flux bins, respectively, in good agreement with 
the CDFs. Based on these number counts of serendipitous 
sources, we find no evidence that the radio sources lie in 



the same cluster en vironments as is observ ed in the case 
of PKS 1138-262 bv lPentericci et alJ <|2002f) . However, as 
indicated in Table the factor ~ 2 excess of faint, > 3cr 
sources in the field of PKS 1138-262 only becomes appar- 
ent in the full 40ks exposure, which may indicate that our 
exposures are not deep enough to make a good compari- 
son. 

To further study whether there is a preference for field 
sources to lie in the vicinity of high-rcdshift radio sources, 
we register the five fields using the radio source posi- 
tions as centroids, and we rotate the fields around their 
original center so that the area of overlap is maximized. 



R. A. Overzier et al.: Chandra observations of z ~ 2 radio galaxies 



15 




Fig. 9. The positions of serendipitous X-ray sources with flux > 1.5 x 10~ 15 erg cm~ 2 s _1 in the five HzRG fields 
combined. Different symbols indicate the sources in the 5 individual fields. The five images are aligned so that the 
radio sources (indicated by large diamond) all coincide with the position of 0156-252, and rotated to maximize the 
overlapping areas of the different pointings (indicated by dotted circles). The full circle indicates an area of 1 Mpc in 
radius around the radio source. 



This is shown in Fig. The number density of sources 
found within 1 Mpc (at z ~ 2) from the radio sources 
is 0.77 arcmin -2 , compared to 0.4 arcmin -2 for the den- 
sity of sources outside this region and 0.6 arcmin -2 for 
the average density of the entire field. However, the cos- 
mic \2™iice_rf field X-ray sources is significant (25%, 
e.g. ICapoi et alJl200lh . making it impossible to conclude 
if serendipitous sources cluster around radio sources on 
the basis of the current data. Spectroscopy may confirm 
whether some of these sources are associated with the ra- 
dio galaxies. 

4. Summary and conclusions 

We have studied X-ray observations of five radio galaxies 
at 2 < z < 2.6, thereby significantly increasing the number 
of Chandra studies on high redshift radio galaxies. The 
main conclusions from our analysis are the following. 

The X-ray emission that we detect from the nuclei are 
consistent with obscured power-law spectra as observed 
for powerful radio galaxies over a wide redshift range (e.g. 
Harris et alJl200d ICarilli et allEool lYoung et all 12002: 
Hardcastle et al.ll2002HScharf et alJl2003h. and can be ex- 
plained by the unified model in which the broad-line region 
of radio galaxies is obscured by a dusty torus surrounding 
the nucleus. For 0156-252, this conclusion is confirmed by 
the fact that the best-fit core spectrum is an absorbed 
power-law with n(HI) ~ 2 x 10 22 cm -2 . 

Extended X-ray emission coincident with the radio 
lobes was detected for several of the sources, albeit at low 
S/N. For all but one source, the straight extrapolation 



from the radio to the X-ray using the radio spectral index 
rules out a synchrotron origin of the emission, unless there 
are large errors in the radio spectral indices measured be- 
tween 5 and 8 Ghz. Although the predicted X-ray syn- 
chrotron flux in source 0156-252 is close to the observed 
value, we have found evidence for the X-ray emission be- 
ing likely associated with shocked, line emitting gas. Our 
observations confirm that the radio lobes of these high red- 
shift sources may interact with the surrounding (forming) 
IGM. 

We interpret the X-ray emission of the remaining 
sources as being due to the inverse Compton scattering of 
CMB photons off radio synchrotron electrons. This con- 
clusion is supported by the fact that our estimates of the 
IC / CMB and equipartition magnetic field strenghts are in 
good agreement. Although IC/CMB is usually dominated 
by other processes (e.g. SSC) at low redshift, the relative 
ease with which it is detected in these HzRGs may be as- 
cribed to the (1 + z) A increase in the energy density of the 
CMB. 

This research was partly motivated by the large RM 
observed for a significant fraction of HzRGs. Can the large 
RM be caused by cluster-sized atmospheres surrounding 
the radio sources? Taking reasonable estimates for the 
pathlength (~ 100 kpc) and the cluster magnetic field 
strength (~ 10/iG), an intrinsic RM of ~ 1000 rad m -2 
would require an ICM electron density of ~ 1 x 10 -3 cm -3 . 
The existence of cores much denser than this is unlikely, 
since they would have an X-ray luminosity of > 10 erg 
s -1 which has not been observed in any of the HzRGs 
observed to this date. Although higher density material 
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(n e ~ 200 cm -3 ) is present in many HzRGs in the form 
of 10 4 K emission line clouds, the filling factor of this gas 
(/ e ~ 10~ 6 ) usually impli ed makes it har d to reproduce 
the high values of RM fe.g. |Pentericcill999(l . Alternatively, 
the large RM can arise from the radio emission passing 
through a sheath of shocked gas surrounding t he radio 
lobes Ce.g. lAthreva et alJll998t ICarilli et aTll2002t and ref- 
erences therein). The increased gas density implied can 
then explain the large RM if the magnetic field is ordered 
on scales of only a few kpc. If the X-ray emission from the 
NE lobe of 0156-252 is from shocked gas (as suggested by 
the bow shock feature seen in the HST/WFPC2 image, 
see Fig. 0), then B ~ 13^G given the RM from Table FJ 
and the derived density of 0.05 cm~ 3 . 

The existence of protoclusters around several HzRGs 
at z r-j 2 and higher has been established, mainly through 
narrow-band, Lya imaging observations and spectroscopy. 
However, such studies have not yet been carried out for the 
sources presented in this paper. An overdensity of X-ray 
sources associated with the well-known galaxy protoclus- 
ter around the radio sourc e 1138-262 at z = 2.2 was found 
bv lPentericci et al] l)2002j) . We have analysed the number 
density of X-ray sources in each of the five Chandra fields. 
None of the fields showed evidence of large-scale structure 
associated with the radio sources. This may suggest that 
these HzRGs are not found in the same, overdense en- 
vironments as 1138-262. Similarly employing several dif- 
ferent methods we found no evidence for virialised gas, 
although the upper limits that we derive on thermal gas 
are not inconsistent with a direct, no-evolution extrap- 
olation of local X-ray luminous clusters out to z ~ 2. 
Based on the X-ray observations presented in this paper 
alone, the current standing is that these five HzRGs are 
not in galaxy overdensities. However, we remark that X- 
ray observations are not the most effective way to search 
for galaxy overdensities given the relatively small number 
fraction of AGN expected at each particular epoch. 

Recent observations of some of the most distant, X- 
ray luminous clusters made with the Advanced Camera 
for Surveys aboard HST, show a color-magnitude rela- 
tion comparable to that of local clusters, indicating that 
early-type galax i es wer e already well-established by z <~ 1 
(Bla keslee etaD Ho03). Thus, the epoch in which the 
global relations that exist in clusters today are shaped 
must probably be sought at significantly higher redshifts 
than currently probed (i.e. z > 1.3). About a dozen high 
redshift galaxy overdensities or protoclusters have now 
been found at z > 2, either around HzRGs and quasars 
or in wide-field surveys. In the near future we may ex- 
pect several candidates suitable for a detailed study with 
Chandra that could determine exactly at what redshift the 
virialised gas was established. However, the sensitivity re- 
quired may be of the order of that of the Chandra deep 
fields, due to the extreme cosmological surface brightness 
dimming. While HzRGs have so far presented the best 
evidence of being associated with massive, forming clus- 
ters, they usually come with a plethora of non-thermal 
X-ray mechanisms. Especially if the radio sources lie at 



the bottom of the potential well that coincides with the 
peak of the bremsstrahlung luminosity, they may not be 
ideal targets for trying to detect this extemely low surface 
brightness emission from the ICM. 

Acknowledgements. The work at SAO was partially supported 
by NASA grant G02-3139B and contract NAS8-39073. We 
thank Michiel Reuland and Andrew Zirm for productive dis- 
cussions, and Melanie Johnston-Hollitt for helping out with the 
MIRIAD software. We are grateful to the referee, P. Tozzi, for 
his many good comments. 

References 

Athreya, R. M., Kapahi, V. K., McCarthy, P. J., & van 

Breugel, W. 1998, A&A, 329, 809 
Bagchi, J., Pislar, V., & Lima Neto, G. B. 1998, MNRAS, 

296, L23 

Band, D. L. & Grindlay, J. E. 1985, ApJ, 298, 128 
Belsole, E., Worrall, D. M., Hardcastle, M. J., Birkinshaw, 

M., & Lawrence, C. R. 2004, MNRAS, 352, 924 
Blakeslee, J. P., Franx, M., Postman, M., et al. 2003, ApJ, 

596, L143 

Brand, K., Rawlings, S., Hill, G. J., et al. 2003, MNRAS, 
344, 283 

Brinkmann, W., Laurent-Muehleisen, S. A., Voges, W., 

et al. 2000, A&A, 356, 445 
Brunctti, G. 2000, Astroparticle Physics, 13, 107 
Cappi, M., Mazzotta, P., Elvis, M., et al. 2001, ApJ, 548, 

624 

Carilli, C. L., Harris, D. E., Pentericci, L., et al. 2002, 
ApJ, 567, 781 

Carilli, C. L., Roettgering, H. J. A., van Ojik, R., Miley, 

G. K., & van Breugel, W. J. M. 1997, ApJS, 109, 1 
Carilli, C. L. & Taylor, G. B. 2002, ARA&A, 40, 319 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Donahue, M., Daly, R. A., & Horner, D. J. 2003, ApJ, 
584, 643 

Eilek, J. A. & Owen, F. N. 2002, ApJ, 567, 202 
Ensslin, T. A., Vogt, C, Clarke, T. E., & Taylor, G. B. 

2003, ApJ, 597, 870 
Fabian, A. C, Crawford, C. S., & Iwasawa, K. 2002, 

MNRAS, 331, L57 
Fabian, A. C, Sanders, J. S., Crawford, C. S., & Ettori, 

S. 2003, MNRAS, 341, 729 
Gehrels, N. 1986, ApJ, 303, 336 

Giacconi, R., Rosati, P., Tozzi, P., et al. 2001, ApJ, 551, 
624 

Hardcastle, M. J., Birkinshaw, M., Cameron, R. A., et al. 

2002, ApJ, 581, 948 
Hardcastle, M. J. & Worrall, D. M. 1999, MNRAS, 309, 

969 

Harris, D. E. & Grindlay, J. E. 1979, MNRAS, 188, 25 
Harris, D. E., Nulsen, P. E. J., Ponman, T. J., et al. 2000, 
ApJ, 530, L81 

Ivison, R. J., Dunlop, J. S., Smail, I., et al. 2000, ApJ, 
542, 27 

Iwamuro, F., Motohara, K., Maihara, T., et al. 2003, ApJ, 
598, 178 



R. A. Overzier et al.: Chandra observations of z ~ 2 radio galaxies 



17 



Kurk, J. 2003, Ph.D. Thesis, Universiteit Leiden 

Kurk, J. D., Rottgcring, H. J. A., Pentericci, L., et al. 

2000, A&A, 358, LI 

Leahy, J. P. 1987, MNRAS, 226, 433 
McCarthy, P. J., Baum, S. A., & Spinrad, H. 1996, ApJS, 
106, 281 

Miley, G. 1980, ARA&A, 18, 165 

Miley, G. K., Overzier, R. A., Tsvctanov, Z. I., et al. 2004, 

Nature, 427, 47 
Mushotzky, R. F., Cowie, L. L., Barger, A. J., & Arnaud, 

K. A. 2000, Nature, 404, 459 
Pentericci, L. 1999, Ph.D. Thesis 

Pentericci, L., Kurk, J. D., Carilli, C. L., et al. 2002, A&A, 
396, 109 

Pentericci, L., Kurk, J. D., Rottgering, H. J. A., et al. 

2000a, A&A, 361, L25 
Pentericci, L., McCarthy, P. J., Rottgering, H. J. A., et al. 

2001, ApJS, 135, 63 

Pentericci, L., Rottgcring, H. J. A., Miley, G. K., et al. 

1999, A&A, 341, 329 
Pentericci, L., Van Reeven, W., Carilli, C. L., Rottgering, 

H. J. A., & Miley, G. K. 2000b, A&AS, 145, 121 
Perley, R. A. & Taylor, G. B. 1991, AJ, 101, 1623 
Rosati, P., Stanford, S. A., Eisenhardt, P. R., et al. 1999, 

AJ, 118, 76 

Rosati, P., Tozzi, P., Ettori, S., et al. 2004, AJ, 127, 230 
Rudnick, L. & Blundell, K. M. 2003, ApJ, 588, 143 
Rush, B., McCarthy, P. J., Athreya, R. M., & Persson, 

S. E. 1997, ApJ, 484, 163 
Scharf, C, Smail, L, Ivison, R., ct al. 2003, ApJ, 596, 105 
Simard-Normandin, M. & Kronberg, P. P. 1980, ApJ, 242, 

74 

Smail, I., Scharf, C. A., Ivison, R. J., ct al. 2003, ApJ, 
599, 86 

Smith, D. A., Wilson, A. S., Arnaud, K. A., Terashima, 

Y., & Young, A. J. 2002, ApJ, 565, 195 
Stanford, S. A., Holden, B., Rosati, P., ct al. 2001, ApJ, 

552, 504 

Stevens, J. A., Ivison, R. J., Dunlop, J. S., et al. 2003, 

Nature, 425, 264 
Taylor, G. B., Barton, E. J., & Ge, J. 1994, AJ, 107, 1942 
Tozzi, P., Rosati, P., Nonino, M., ct al. 2001, ApJ, 562, 42 
van Ojik, R., Roettgering, H. J. A., Miley, G. K., & 

Hunstcad, R. W. 1997, A&A, 317, 358 
Venemans, B. P., Kurk, J. D., Miley, G. K., et al. 2002, 

ApJ, 569, Lll 

Wilson, A. S., Young, A. J., & Shopbell, P. L. 2000, ApJ, 
544, L27 

Windhorst, R. A., Keel, W. C, & Pascarelle, S. M. 1998, 
ApJ, 494, L27 

Young, A. J., Wilson, A. S., Terashima, Y., Arnaud, K. A., 

& Smith, D. A. 2002, ApJ, 564, 176 
Yuan, W., Fabian, A. C, Celotti, A., & Jonkcr, P. G. 2003, 

MNRAS, 346, L7 



